Tension-induced non-linearities of fiexural modes in nanomechanical resonators 
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We consider the tension-induced non-linearities of mechanical resonators, and derive the Hamil- 
tonian of the fiexural modes up to the fourth order in the position operators. This tension can be 
controlled by a nearby gate voltage. We focus on systems which allow large deformations u(x) 3> h 
compared to the thickness h of the resonator and show that in this case the third-order coupling can 
become non-zero due to the induced dc deformation and offers the possibility to realize equations of 
motion encountered in optomechanics. The fourth-order coupling is relevant especially for relatively 
low voltages. It can be detected by accessing the Duffing regime, and by measuring frequency shifts 
due to mode-mode coupling. 

PACS numbers: 85.85. +j 



Recent progress in fabricating nanomechanical res- 
onators has shown how these systems can be used for 
ultrasensitive measurements of mass, force and charge 
PHI]. Within the couple of past years these systems 
have also entered the quantum realm [5] as superposi- 
tions of vibration states and zero-point vibrations have 
been measured. Even though such measurements can be 
performed in a regime where the elastic properties of the 
resonators could essentially be considered as linear, the 
extension to non-linear conditions is well within reach of 
the current experimental techniques. 

In this paper we consider the generic non-linearities 
of the resonators, how these show up in measurements, 
and how they arise when the resonators are manipulated 
electronically. In general, the effect of non-linearites 
is twofold: on one hand they modify in an amplitude- 
dependent way the resonant frequency of a given normal 
mode (Duffing self-non-linearity); on the other, they in- 
troduce a coupling between normal modes. Such non- 
linearities show up in the presence of strong external 
driving, which allows to control the coupling of different 
modes or to detect their occupation numbers. 

Motivated by the recent advances in fabricating 
graphene and carbon nanotube resonators [H [3] , we con- 
centrate especially on the regime of thin resonators where 
the mechanical deformation can be large compared to the 
resonator thickness. In this case, the major source of non- 
linearity is the tension induced by the deformation itself. 
Starting from the mechanical energy of the deformations, 
we derive the generic Hamiltonian of the fiexural modes, 
including non-linearities up to the fourth order in the 
vibration amplitudes. We assume that the dc deforma- 
tion, dictating the strength of non-linearity, is driven by 
a nearby gate voltage as in Fig. [T] Concentrating first 
on the Duffing self-non-linearity of the modes, this then 
allows us to derive the voltage dependence of the Duffing 
constant and show that it changes sign for a certain value 
of voltage that depends on mode index and the amount of 
initial tension. Therefore, studies of the Duffing constant 
reveal information about the parameters of the system, 
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FIG. 1. Schematic picture of the studied setup : a metallic 
beam whose deformation is controlled by a gate voltage V g 
coupled to the beam via capacitance C(d — u(x, t)). 

in particular on the initial tension, which may otherwise 
be difficult to obtain by only concentrating on the voltage 
dependence of the mode eigenfrequencies. We go on to 
analyze the inter-mode coupling and show that the non- 
linearities allow creating a radiation-pressure-type cou- 
pling between the different fiexural modes. Such a cou- 
pling allows realizing optomcchanics-type experiments, 
where one of the modes is cooled or heated by driving 
another mode. We provide quantitative predictions for 
the optical spring effect (driving-induced frequency shift) 
and changes in effective mode damping responsible for 
the cooling/heating behavior and show how these can be 
tuned by the dc gate voltage. The general Hamiltonian 
describing a non-linear resonator is of the form 
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Here x r , 



a n are the dimensionless position op- 



erators. The non-linearities are described by the coeffi- 



cients TV, 



and F~, 



ilk • 



In the following we consider a 



mechanical resonator exhibiting the non-linearities dis- 
cussed above. We analyse a beam of mass to with length 
L, thickness h and cross-section S suspended on top of 
a gate capacitor at voltage V g (see Fig. [T]) . The flexural 
vibrations are characterized by the deformation u(x, t) of 
the beam. Defining z = y/L, u! = d z u and introducing 
the notation (u\v) = f Q u(z,t) ■ v(z,t) dz one can obtain 
|TJ from the elastic energy of a resonator [7] 

e[u(x,t)]=^(u"\u")+ g («>') 



ES 

+ ^ (u'W) 2 + s {oIce {u(x,t)], (2) 



with E the Young modulus, I y the bending moment, To 
the initial tension of the resonator and £f or ce the poten- 
tial energy of the force acting on the resonator. The 
latter is of the form £f orco = — (V 2 /2) f Q C[d — u(z,t)]dz, 
where C[d] is the capacitance between the gate and the 
beam at the distance d. In order to arrive at a Hamil- 
tonian of the form given in Eq. 0, we assume that 
the gate voltage V g is the sum of a DC part Vdc and 
a small AC part V ac -C Vdc- These voltages lead to a 
static deformation Ud c (x) <C d and a time varying part 
u ac {x,t) <C Udc(x). Expanding u ac (x,t) on an arbitrary 
basis Xn(x), u ac = h J2 n V n (t)Xn(x) , w e can write the po- 
tential energy containing terms with two, three and four 
y n s. Writing the Hamiltonian in terms of the stress en- 
ergy Stress = muj 2 s h 2 /2 = ESh 4 /(8L 3 ), these terms are 
characterized by the dimensionless parameters which we 
denote by [fi 2 ]^, A™ lo , and 6™^ for the second-, third- 
and fourth-order terms, respectively. These parameters 
are described in detail in the Appendix [8 . In particular, 
they depend on the dc bias voltage Vd c and the total ten- 
sion T in the beam. As discussed below, the latter also 
depends on Vdc- The behavior of the coefficient [£1 2 ]™> 
determines the voltage dependence of the eigenfrequency 
as described in [9l[T0]. For Vdc = 0, the third-order terms 
AJ^ vanish because of symmetry but for large Vdc they 
grow as AJ^ D cx V 2 / 3 . The voltage dependence of the 
fourth-order terms O no on the other hand is weak and in 
our analytical approximations [5] disregarded altogether. 

We arrive at the desired form by writing the Hamil- 
tonian in a basis which diagonalises [^ 2 ] m and scaling 
the amplitude hy n of each mode by its zero point motion 

-i— [llj . Non-linearities are coming primarily 
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from the induced tension, which is maximized for large 
deformations u(x) 3> h. Therefore we concentrate on sys- 
tems which allow large deformations, i.e., systems with 
c(> h. We truncate the expansion to the fourth order as 
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FIG. 2. (Color online) Frequency response function for the 
first mode with 7 2 = Q.Qlui'i and \f d \ 2 = 0.002 m 2 uj 2 s h 2 . Here 
-Bstross = ESh 4 1 '(8L 3 ). 



in the model employed above the higher-order couplings 
are relevant only close to the point where the beam pulls 
into contact with the gate plane [TU] , 

As an example, a single-layer graphene sheet with h = 
0.34 nm, L — 1 /im, S/h = 1 /im, E = 1 TPa and mass 
density p — 1400 kg/m 3 would have to s = 250 MHz and 
E stICSS = 0.01 eV. 

Self-non-linearity. Let us first consider the non-linear 
effects which occur when driving mode n with a driving 
force fd cos(ojrft), uid ~ w„, disregarding the coupling to 
the other modes, as in the absence of direct driving of the 
other modes these would show up only in a higher order 
in the non-linear coupling constants. We also exclude 
the special case when 2 or 3 times the mode frequency 
does not match one of the other mode frequencies [T2T 
I14j . Including dissipation the equation of motion for the 
amplitude of mode n is 

i n + w 2 x„ + jx n + 3T n x 2 + 4J r „.T 3 1 = — cos(w d t). (3) 
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T n = hjfi (0™") d and the subscript d denotes that the 
tensors are written in the basis which diagonalises [f2 2 ]. 
The frequency response equation of ^ can be solved 

from [rjnn] 
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where 7 is the damping constant and D 
f(k) 2 = ^ 2 /(2/ l 2 )[(e««) d -5(A«J 2 /[2(^) 2 ]]. 
This frequency response function is the same as one 
would get when considering only the fourth-order non- 
linearity T n , i.e, a Duffing oscillator. The effect of the 
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third-order non-linearity T n is to shift the value of the 
cubic non-linearity in the frequency response function. 

An example response function obtained for different 
dc gate voltages is plotted in Fig. [2] and shows up a 
crossover from D > (hardening) to D < (soften- 
ing). The dc voltage dependence of D is plotted in the 
Appendix [5] for a few different magnitudes of initial 
tension To, characterized by the dimensionless quantity 
t = AT L 2 /(ESh 2 ) = h 2 T /(2E strcss L). Here we char- 
acterize its overall behaviour. The behavior of the Duff- 
ing constant depends on the total tension T of the beam 
which is the sum of the initial tension To and the ten- 
sion induced by the deformation Ud c {x) caused by Vdc- 
The latter has to be calculated self-consistently from the 
Euler-Bernoulli equation as discussed in [9]. In what 
follows, we describe this behavior in terms of the dimen- 
sionless quantity r = ATL 2 /(ESh 2 ) = h 2 T/(2E strcss L). 
In the limit d ^S> h it satisfies [5] 
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(— hddCV 2 ) / E stress . This equation is valid 
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provided the resultant t > 1. In the same limit, we find 
for the Duffing coefficient of the fundamental mode n = 1 



D = 



2h 2 



2 V 
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where /„(r) = 2nvr[2/(n7r) 2 - V3r/((n7r) 2 + 3r)] for 
n odd and zero otherwise. For a rectangular beam 
I y /{Sh 2 ) = 1/12, but the overall behavior of D does 
not greatly depend on the exact shape of the beam. 

Solving Eqs. ^ and Q allows us to find the approxi- 
mate behavior of the Duffing constant as the gate voltage 
is tuned. We find that D changes sign at a value of the 
gate voltage V* that can be quite well fitted to the func- 
tion V* « %/2r 3 / 4 + 8 (see 0) or 
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Moreover, at large values of the voltage V g ^> V* , D (ap- 
proximatively) saturates to the value D sat ps — I20u 2 /h 2 . 

Contrary to the fundamental mode n — 1, the deforma- 
tion induced changes in the Duffing constant of higher- 
order modes are rather small compared to its value for 
V g = 0. 

Non-linear mode coupling. Let us now concentrate on 
the non-linear coupling between the modes [16]. The 
regime investigated here is formally analogous to the 
setup encountered in optomechanical systems [T7^l20j . 
where an external driving electromagnetic field, coupled 



to a resonant cavity, alters the characteristic response pa- 
rameters of a mechanical resonator. More specifically, by 
aptly tuning the pump frequency, it is possible to alter 
the resonant frequency of the mechanical resonator (opti- 
cal spring effect [2"UH2"2"] ) and its damping, thereby induc- 
ing cooling [T!5] or amplification |20j . Here we consider 
the case where one mechanical mode, say with eigenfre- 
quency u mi corresponds to the cavity mode, and another 
one, uj n , to the mechanical mode. We also assume that 
uj m > 2uj n . Let us discuss what happens if the system 
is driven with a frequency Ud — ui m — A, A ±w„ and 
probed around uj n . 

Neglecting other modes, the Hamiltonian is of the form 

H = u n a\d n + bj m a\a m + T n x\ + T m x 3 m 

T T nnrn x n x m -4- T nmm x n x m -4- F n x n -\- F m x rn 
~t~ F nnnrn x n x m -\- F nnmrn x n x m T F nmrnnl x n x rn: (8) 

where T nmo and F nmop are the sum of all the permuta- 
tions of indices n, m, o, p of T„ mo and ¥ nmop , respec- 

tively, and T nmo — r x Z p n x Z p m x Z p o (A mo )^, IF nmo / — 

Tinn and T m = 



,, x zp m x zpi x zp„ (®mi)d> ^™ 
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Tmmm- Using the input /output formalism [23] the equa- 
tions of motion for operators a n and 



u n d n + 4F n xl + 3T n xl 



(9) 
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T 2(T mmn x n -\- F nn7nm x n ^d n 
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Here we have written operator a m in a frame rotating 



with frequency ujd- We linearize (10) and (11), rewriting 



the operators as a sum of a static a and a fluctuating part 
(5a, a = a + Sd. Keeping terms which are of the order of 

l > am > . Solving for 



-4Tn 



(^ a -) we obtain a n + a* n 
5a n we find the frequency response function for the input 
signal <5a™. It is a Lorentzian function peaked at 



4|a m | 2 



t t i r 

p r. ±n-Lnmm, ^ ± nmm 



2 w„ 



and whose width is 

7n,eff — 7n 



4a r 



±- 



T 2 



Qr, 



(12) 



Here Q m = uj m /"f m is the quality factor of mode m and 
we have assumed for simplicity the fully side-band re- 
solved limit uj n 3> 7 m . We remark that the results are 
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similar to those obtained in optomechanics, the only dif- 
ference comes from the second term in the effective fre- 
quency which is proportional to the self-non-linearity. As 
in the case of Duffing non-linearity we consider the limit 
of r > 1. We find that the effective frequency when driv- 
ing mode m = 3 and probing mode n = 1 depends on 
the gate voltage as 



Wl.eff — LO\ _ 

!)\~ ; = 
where g n (V,r) 



3/i(r) 



2V" 
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2 V 4 
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2 V" 



(13) 



and 



2„2 

ZPn 



-p 1 describes the amplitude of the 
pump. The effective damping changes as 



7i,cff-7i =± 2(3tt) 4 / 1 (t) 
xQa ~ 93(V,t) 



4 
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In Fig. [3] we plot the effective frequency and the effec- 
tive damping when driving mode m = 3 and probing its 
effect on mode n = 1. We see that both for the red- 
and blue-detuned pumping, i.e., A = ±wi, the frequency 
shift induced by pumping, uji^r — lj-i .[23] is positive at 
low gate voltages due to the fourth-order term Fn 33 in 
Eq. (11), changes sign upon an increasing dc gate volt- 
age, and tends to a voltage-independent value at large 
voltages. The fact that the overall frequency shift is in 
both cases negative — in contrast to the traditional op- 



tomechanics — results from the second term in Eq. (11 1, 
which reflects the effect of the self-non-linearity T n , and 
which is independent of the sign of A. This behavior ap- 
plies only to the combination n = 1, m = 3. For higher- 
order n, the voltage-induced changes are small compared 
to the frequency shift at Vd c = 0. However, choosing 
n = 1 and higher m results into more complex behavior 
and the spring effect wi. c fr — u>i may change sign more 
than once in the case of blue detuning (see [3]). 

On the other hand, the change in the effective damping 
(inset of Fig. [3]) depends on the sign of A. For red detun- 
ing, A = w n , 7„ i0 ff increases as the voltage is increased, 
whereas for blue detuning 7„ ie ff decreases. For a fixed 
amount of fluctuations coupling to mode n, the increase 
in damping leads to (side-band) cooling [TS] , whereas the 
decreasing damping leads to heating and, when j e g be- 
comes zero, to a parametric instability 25J. Between 
these regimes, the blue-detuned driving can be used for 
signal amplification |20) . 

In conclusion we have derived the Hamiltonian of a 
thin doubly clamped nanomcchanical resonator taking 
into account the non-linearities between the amplitudes 
of the flexural modes induced by a nearby gate voltage. 
Besides the Duffing non-linearity, we also find a third 
order non-linearity directly related to the DC deforma- 
tion of the beam. This third order non-linearity adds 
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FIG. 3. (Color online) Effective frequency and damping (in- 
set) of mode n — 1 when driving mode m = 3 with initial ten- 
sion ro = (no symbols) and r = 10 (circles) in the case of red 
detuning (red lines, lower) with A = uji and in the case of blue 
detuning (blue lines, upper) with A = — The full lines are 
numerical results obtained by solving the full Euler-Bernoulli 
equation obtained by requiring u(x, t) to minimize the energy 
in and dashed lines follows Eqs. |i"3] ) and (25 I. Here 



h 2 



and E s - 



ESh A /(fiL 3 



to the Duffing non-linearity and changes the behavior 
of the frequency response function. Besides the self-non- 
linearity described by the Duffing behaviour, we find that 
the different modes of the beam are non-linearly coupled. 
The effective Hamiltonian of a pair of such modes resem- 
bles that of a mechanical degree of freedom coupled to a 
cavity, with the difference that in the current setup the 
cavity is replaced by another flexural mode. Therefore, 
such a coupling offers the possibility of observing the mo- 
tion of one mode by observing its effect on another mode. 
Such effects are the spring effect and the changing damp- 
ing, and the latter can be used for side-band cooling or 
amplification of a given mechanical mode. 

Besides the non-linearity induced by bending described 
here, there may be other sources of non-linearity in thin 
metallic beams, such as those related with non-linearities 
in electronic properties [2 6) or non-linearities induced 
by stretching. Our results help to identify the direct 
bending-induced non-linearities and therefore facilitate 
of a precise tuning of nanomcchanical resonances. 
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APPENDIX : DERIVATION OF THE NON-LINEAR HAMILTONIAN 

Expanding u ac in Eq. (2) of the main text on an arbitrary basis Xn( x )j u ac = n DnXn we get the Hamiltonian 
of the form 



E — ~t~ -Estrcss 

2m 



{ [ fi2 ] 1 Vny m + KnoVny m y° + QZ P yny m y y p } + E g , 



(15) 



with the non-linear coefficients 



[ fi2 C = - 2 (X'L\X'L) +rUnWJ +2(x'„\x' m )«M dc ) +^u' dc \x' n ){u' dc \x' m 



A" 



QUO 

w nip 



dc 



2 h 2 



1 Al 



d 2 C\d- 



i(x,t)) 



du 2 r . 



^' dc \x' ){x'n\x' m 
(Xn\x'm)(X o I X.p) 



vl 



{Xn\Xm)dx 



1 d 3 C[d~ u(x,t)] 



muj 2 h 2 J dul c 

1 d 4 C[d~ u(x,t)} 



XnXmXodx 



Vr 



mu> 2 h 2 J 



dut 



XnXmXpXodx 



E g = (2V dc V ac + V 2 C ) / f(u (x))dx. 
Jo 



(16a) 
(16b) 
(16c) 
(16d) 



Here C[d] is the capacitance between the gate and the The solution of this integro-differential equation is [9] 
beam at the distance d, \ 



\h 2 = 



ergy of a beam displaced by h and 



EI 
2L3 



is the bending en- 

2 h 2 = ES£ ig the 



2 ""*'«'« ~ 8L 3 

stress energy of the beam displaced by h with respect to 
its equilibrium position. The coefficient E g describes the 
feedback of the motion of the resonator on the gate volt- 
age and is neglected below as we assume a fixed voltage 
drive. Besides the voltage, the system is described by the 
two dimensionless parameters to = ATqL 2 / (ESh 2 ) and 
(ljo/u s ) = 4I y /(Sh 2 ). For a rectangular beam, which 
we consider in the following, (uj /u! s )' 2 = 1/3. Overall, 
our main results do not greatly depend on wo- The thick- 
ness h appears in the above expressions only because it 
sets the magnitude of the deformation — it scales out 
from the final results of observable quantities. 

Writing the Hamiltonian in the basis which diago- 
nalises Q and scaling hy n by the amplitude of the zero- 
point motion \J 1 /mu) n one arrives at Eq. (1) of the main 
text. We consider here specifically a doubly clamped 
beam. For low voltage the DC deformation u dc is given 
by the Euler-Bernoulli equation in the case of a parallel 
plate capacitance model: 



mui Q hu dc — [ moJ T h + 2mui s h 



'dc: 



dx 



V 2 eWL 

~2 dP~' 
(17) 



v 2 ( (r, 

u dc = I coth I - ) (cosh(£:r) - 1) - sinh(£x) 



(18) 



+ £x - , 

where V 2 = (-hd d CV 2 )/E stTCSS , E stTess = muj 2 h 2 /2 
ESh 4 /(clL 3 ). We define the total tension 



t = t + 2 / u dc dx 
'o 



(19) 



and £ = V^t- Substituting Eq. ^ into Eq. (TTsl) , 



integrating and then disregarding exponential terms 
exp(— £) yields a self-consistency equation for r, 



T = T 



96r 2 



3V3 



'T T 



(20) 



Since in the limit of large t Eq. (17) reduces to the wave equation, we use the harmonic wave function Xr. 
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\/2 sin(n7ra;) as a basis for (17) and get 
(rnr) 4 



\n 2 ] n 

L J m 

A" 



QUO 



+ T(^) 2 + ^/„(T)/ m (r) 



~2~d 



V 2 



V2(nn) 2 — f (T)5Z 



V 2 



2r' 



2 

2 / l \ 3 



4^7^710 ((-1)™+™+° - 1) 



7r(m — n — o)(m + n — o)(m — n + o)(m + n + o) 

2 2 ri o V r2 //lVl 
(7171") (j57r) <5 m <5 p — l^fin+m+o^-p ^n+T»,a+P + ^n+p,m+o + ^n+o.m+p — $n 



+m+o,p 



&n+m-\-p,o ^n+o+p,m ^m+o+p, n\ 



/ n (r)=n7r(l-(-lD 



'3t 



7r 2 n 2 7r 2 n 2 + 3t 



(21a) 
(21b) 

(21c) 
(21d) 



The last terms in Eqs. (21a), (21b) and (21c) are relevant 



only for low voltages in the presence of initial tension, and 
do not greatly contribute to the physics discussed in this 
paper. We thus drop them out in the ensuing analyt- 
ical approximations. However in the numerical results, 
we use the full solutions of the Euler-Bernoulli equation 
to determine the eigenmodes and the coupling constants. 
Nevertheless, Eqs. ( 21a||21c | represent a fairly approxi- 
mations in the limit of relatively strong tension. 
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FIG. 4. Duffing constant for the first mode with different 
initial tension r. From left to right the dimensionless initial 
tension to starts from and increases with the step of 10. 
The dashed lines are our analytic expressions and the full line 
are numerical solutions of the full Euler-Bernoulli equations 
obtained from Eq. (2) of the main text. The deviation be- 
tween the two set of curves at low V is due to our scheme of 
approximating mode functions by harmonic functions. 
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FIG. 5. Duffing constant for the third mode with different 
initial tensions to. From left to right the dimensionless initial 
tension tq starts from and increases with the step of 10. 
The dashed lines are our analytic expressions and the full line 
are numerical solutions of the full Euler-Bernoulli equations 
obtained from Eq. (2) of the main text. The approximation 
with harmonic mode functions underestimates 0™ p , which 
is the reason for the discrepancy between the full numerical 
solutions and our analytic approximations. 



disregarding the off-diagonal terms in Eqs. ( 21a|21b ) and 
solving the tension r with Eq. (20). The tension r 



ex- 
hibits a rather complicated voltage dependence, however 
its behaviour can be investigated in different limiting 
cases. A characteristic value for the voltage can be found 
by substituting t = tq into Eq. ( 20 1 and comparing To to 
the second terms of the right hand side of Eq. ( 20 1 . This 
yields 



V* = 96 1/4 - 



3/4 



1 - 3 



x/3 - 



1/4- 



(22) 



We can estimate the overall behavior of the Duffing for V ^> V* we have t 



Thus we find that for V -C V* the tension r 



constant D = uj 2 /(2h 2 ) [(9 



nn \ 
nn)d 



96 1 / 3 1 



t q while 
Disregarding terms 



(A™„) 2 /[2(S!™) 2 ]J by coming from the electrostatic force the general expression 
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FIG. 6. Crossover voltage V c for the sign change of the Duffing 
constant with respect to the initial tension r. 



for the Duffing constant is 
5 



D 



2h 2 



(wr) 4 



(«7r) 4 / n (r 



,2V 1 



8§K)+rW + / n (r)^ 



(23) 

As shown in Figs.|4j|5j at low V, D starts from a positive 
value D(V = 0) « u; 2 / (2h) 2 (nir) 4 and tends to a voltage- 
independent value D(t0 = u 2 J(2h 2 ) [ ( ""f 9 [ 2 ( ^£~ 768) 

at large voltages. Note that, for a symmetric dc defor- 
mation, this behavior is only valid for odd-order modes 
(with symmetric eigenf unctions with respect to the cen- 
ter of the beam). Indeed, from Eq. (21d), we find that 
the Duffing constant is voltage independent for an even 
n. In Fig. |4j we plot the behaviour of the Duffing con- 
stant for the first mode for different values of r. We find 
that D changes its sign for a given value of the voltage, 
V = V c and the effect of the initial tension is to shift 
the crossover voltage to higher values (Fig. [6]). For n > 1 
the Duffing constant varies only weakly with V and stays 
always positive (see Fig. [5j 



MODE COUPLING 

The frequency response function / for the input signal 
5a™ solved from Eqs. (10-11) of the main text is 



(24) 





\ 2 + {l f - luJ ) 2 -c 2 \a m \ i 


[a + + lu + uj n ] - 2b 2 \a m \ 2 (A + c\a rn \ 2 ) 


A 2 + (2 f - lw ) 2 - c 2|a m | 4 " 




(3f- — iw) 2 + cu 2 + 2aw n 


-4 Wl 6 2 |a m | 2 (A + c|a m |2) 



with 



b 4^ nnmm (o; n -f- cx n ) -\- c 2T nrarrL 
c = \2F m . 

This frequency response function describes a Lorentzian resonance around an effective frequency tv nie ff with damping 



7n,eff described in Eqs. (12) and (13) of the main text. Using the approximations leading to Eqs. (16 1 we find the 
effective frequency 



w n,eff w n 2 2 4 
' = U m 7T 

X 



An 2 



] /(^+T(nn) 2 + fn (T) 2 £)( 



™I+ T ( m7r) 2 + /m(T )2^ 



and effective damping 



7n,cff — In 



= ±- 



2{mnYf n (rf 



-It 2 



xQ m lM^ + r(m7r) 2 + /m(r)2 |^ 



(25) 



Here x = 4mw 2 |a m | 2 a; 2 Pm -jfjf- describes the amplitude of the pump. 
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In Fig. 3 of the main text, we plot the effective frequency and the effective damping when driving mode m = 3 and 
probing its effect on mode n = 1. We see that, for red-detuned pumping with A = w n , the frequency shift induced by 
pumping is positive at low gate voltages, u! n . c s(V g = 0) — oj n = (nir) 2 (mir) 2 x, changes sign with increasing dc gate 
voltage, and tends to a voltage-independent value at large voltages 



w n,cff 



V>Tn,l , , 2 4 
n > X[nm) 7T 



4608/(n7r) 



384/(n7r) 2 



(nvr) 2 + 192/(mr 2 ) n 2 y/[(nn) 2 + 192/(mr 2 )] [(m^) 2 + 192/(tott 2 )] 



(26) 



On the other hand the effective damping increases with an increasing gate voltage until it reaches a voltage- independent 
value 



ln,eS - 

xQr. 



7 _m° 768?r 4 

7r 4 m 4 + 192' 



For blue-detuned driving, when A = — u n , the optical spring effect (w nje fj 
voltages to the value 



W„ nff 



■U3 n V X( nm ) 7T 



4608/ (nit) 



(27) 

uj) increases until it saturates at large 
384/(mr) 2 



(titt) 2 + 192/(mr 2 ) n 2 ^[(nn) 2 + I92/(n7r 2 )] [(tott) 2 + 192/(m7r 2 )] 



(28) 
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FIG. 7. Effective frequency and damping (inset) of mode n = 
1 when driving mode m = 5 and when initial tension to — in 
the case of red detuning (red lines, lower) with A = wi and in 
the case of blue detuning (blue lines, upper) with A = 
The full lines are numerical results obtained by solving the 
full Euler-Bernoulli equation and dashed lines are analytical 
results derived in the text. 
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FIG. 8. Effective frequency and damping (inset) of mode n = 
3 when driving mode m = 5 and when initial tension to — in 
the case of red detuning (red lines, lower) with A = ljs and in 
the case of blue detuning (blue lines, upper) with A = — u>3. 
The full lines are numerical results obtained by solving the 
full Euler-Bernoulli equation and dashed lines are analytical 
results are analytical results derived in the text. 



while the damping decreases until it reaches the value 



7n,cff - 7 
XQm 



768tt 4 



n 4 7T 4 m 4 + 192 



(29) 



These predictions are compared to the full numerical so- 
lutions obtained from the Hamiltonian Eq. (2) of the 
main text in Fig. 3 of the main text. Although we are 
focusing only on the first and third mode in the main 
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text one can also consider pairs of other modes. We plot 
the effective frequency and effective damping for n = 1 
and m = 5 in Fig. [7] and for n = 3 and m — 5 in Fig. [8] 
In the first case when driving mode m = 5 and measur- 
ing mode n — 1 the change in the effective frequency is 
larger than the one we have when driving mode m = 3. 
We also find that in the case of blue detuning the effec- 



tive frequency changes its sign twice as a function of the 
gate voltage. When driving mode m — 5 and measuring 
mode n = 3, the voltage dependence of the change in 
the effective frequency is relatively small as the fourth- 
order term dominates throughout the interesting regime 
of voltages. 



